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Abstract 

One of the most important issues of quantum engineering is the construction of low-dimensional 
structures possessing desirable properties. For example, in different areas of possible applications 
of the structures containing quantum wells (QW) , there is need to have QW energy spectrum, 
which is predetermined. Then the following question arises: can one reconstruct the shape of 
QW which supports this spectrum? We outline the possible strategy of the QW potential shape 
reconstruction, if the spectrum of QW is given in advance. The proposed approach is based on 
the combination of different techniques such as Inverse Scattering Problem Method, Darboux and 
Liouville transformation. It enables to take into account the space-variable dependent effective 
mass of charge carriers and allows the kinetic energy operator to be of non-Hermitian as well 
as Hermitian form. 

The proposed technique allows to construct phase-equivalent potentials, to add the new 
bounded states to (or remove some of them from) the spectrum supported by an initial poten- 
tial and provides a systematic procedure for generating new exactly solvable models. 

PACS:03.65.-w, 03.65.Fd 
Key words: Quantum well, Inverse Scattering Problem, intertwining operator technique, Dar- 
boux transformation, Liouville transformation. 

1 Introduction 

In recent years we have been witnesses of the rapid progress in nanoelectronics which is already 
on the way to continuing the outstanding successes of microelectronics. This became possible 
among others, due to the development of technologies and techniques, such as Molecular Beam 
Epitaxy (MBE) for instance, which enables to deposit thin layers of different materials one on top 
of the other, with almost atomic precision. The last one, in its turn, is able to produce a variety 
of low- dimensional structures, ranging from a heteroj unction formed at a single interface, through 
quantum wells (QW) to superlattices. It would not be an overstatement to say that a new paradigm 
of electronics emerged, for which even the name has been already coined, Quantum Technology or 
Quantum Engineering. It seems, however, that quantum engineering in its present stage, in spite 
of all its successes and maturity, is still passive in the sense that it makes use of, figuratively 
speaking, less 'degrees of freedom' than it possibly could. It means that the 'palette' of QW- 
potential shapes is still limited to a few most popular ones: rectangular, parabolic or semi-parabolic 
and this circumstance obviously restricts the possibility to choose and control the energy spectrum 
of QW produced by means of MBE. Meantime, in different areas of possible applications of the low- 
dimensional structures mentioned above, there is often need to have some specific kind of spectrum 
known beforehand and a question arises: how to produce the QW with a predetermined spectrum? 
At least from the theorist's point of view, this question can be reformulated as follows: suppose 
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the spectrum of QW is known; that is, one knows the number of quantum levels and the distances 
between them (remember, the spectrum is not necessarily equidistant as in case of a parabolic 
well or similar to the spectrum of a square well); can one reconstruct the QW-potential which 
supports this spectrum? An affirmative answer to the question would make quantum engineering 
more flexible and active, providing the opportunity to develop multitude of novel quantum devices. 

The aim of this paper is to develop an approach to the QW-potential reconstruction which 
combines different techniques, such as Inverse Scattering Problem Method, Darboux (or intertwining 
operator technique) and Liouville transformations in reference to the generalized (it means, with 
space-variable dependent mass) Schrodinger equation. 

2 Necessary preliminaries 

As we already mentioned, one often needs to know, how to reconstruct the potential provided that 
the spectrum supported by this potential is known beforehand. One can look at the problem also 
in this way: suppose we have a set of "scattering data", whatever it means; we do not define them 
precisely right now. Can we find the potential which produces this set of data ? Affirmative answer 
to this question was obtained for the first time by mathematicians I.M. Gel'fand, B.M. Levitan 
[1,2] and V.A. Marchenko [3] (see also Refs. therein). They elaborated the method of potential 
reconstruction by means of spectral or scattering data which is now known under the name of 
Inverse Scattering Problem Method (ISP). We shall term their method or approach also as a GLM- 
method (approach). In a particular case of QW, the question is reduced to the following: suppose 
we have a number of quantum levels with the given distances between them; what the potential 
shape (or which the form of) QW should be in this case; or in other words, can we find the potential 
of QW, if the spectrum and the depth of QW are given? In Ref.[4] the provisional positive answer 
to this question was obtained by means of the ISP-method. The reader also can find in Ref. [4] the 
details concerning the motivation for such inquiry. Here, for the sake of consistency and reader's 
convenience, we outline the main ideas of GLM-approach only very briefly. 

Consider the eigenvalue problem of the Schrodinger equation with the potential V(x) in one 
space dimension 

[-^ + V{x^{x,k)=k 2 4>{x,k), (1) 

(here we assume h 2 /2m = 1). The GLM-method may be viewed as a dispersion theory for the 
Schrodinger wave function of (1). ^From the solution to (1) with k complex, one can define Jost 
functions ff and f% which are analytic in the upper-half fc-plane with the following asymptotic 
behaviour 

ff{x,k) ~ exp(±ikx), as x — > +oo, 
f^(x,k) ~ ex.p(^ikx), as x — > — oo, 
and construct a meromorphic function $(x,/c) as 

<E>(x, k) = a _1 (fc)/^"(x, k) exp(ifcx), Im/c > 0; 

= fi(x, k*) exp(ikx), Imk < 0, 

where a~ 1 (k) is the conventional transmission coefficient. 

Thus, <E>(x, k) is completely determined by its singularity structure which consists of a cut along 
the real fc-axis and some number, say N of bound-state poles on the positive imaginary axis. The 
spectral weight of the cut is essentially a scattering-state wave function multiplied by the reflection 
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coefficient, both evaluated at real k. Similarly, the pole residues are essentially constants times 
bound-state wave functions. Upon Fourier transformation the dispersion relation for $ becomes 
the Marchenko integral equation [3] (see also [5] for details) which determines the wave functions. 
The inverse problem of the potential reconstruction by means of scattering data can be reduced, 
generally speaking, to solving an integral equation 



which is given entirely in terms of the reflection coefficient and 2N bound-state parameters. Here 



where S(k), E n and M% are scattering data: S(k) is a scattering matrix, E n = (in n ) 2 are bound- 
state energies and M 2 are its normalization constants. 

If the reflection coefficient can be represented by rational functions of k, the Marchenko equation 
can be solved exactly by algebraic technique. The most simple case corresponds to a reflection 
coefficient vanishing for all real k (S(k) = 1); then the expression for Q(x, x') contains only the sum 
over the bound states, and the integral equation reduces to a system of N linear algebraic equations 
with respect to the kernel K(x,x'). The potential V(x) also can be reconstructed by means of 2N 
parameters. The half of the parameters makes N bound-state energies E n , n = 1,2,...N, while 
the others are the normalizing constants M 2 which in Nuclear Physics are assumed to be obtained 
from scattering data as M 2 = iResS(k)\k = in n . It corresponds to the transparent and symmetric 
potentials. These 2N numbers supply a complete set of the parameters for solving the inverse 
problem to exist and be unique. Obviously, in case of QWs the second half of the scattering 
data (the normalizing constants) is inaccessible. It turns out however, that the way out actually 
exists. The possibility to reconstruct the infinitely deep and symmetric potentials solely by means 
of bounded states was established for the first time by B.M Levitan and M.G. Gasymov who 
proved the The Theorem about Two-spectrum [6,2]. The meaning of this theorem in brief is the 
following: to reconstruct symmetric potentials V(x) = V(—x), it suffices to know only the positions 
of energy levels without knowing the normalizing constants. Later on, this fact was re-discovered by 
H.B.Thacker, C.Quigg and J.L.Rosner [7] who modelled the confining potentials binding massive 
quarks and antiquarks in meson systems. They were able even to fit approximately the lower parts 
of infinitely deep wells (parabolic, harmonic oscillator, linear and rectangular wells) by means of 
a limited number of levels whose positions are known, using some additional technical trick which 
is termed a stabilizing parameter. This observation was used in Ref.[4] in order to reconstruct the 
shape of QW provided that its energy spectrum and depth are known in advance. The problem 
however is, that in the approach developed there, the electron effective mass was supposed to be 
constant and independent on the space variable, while in reality it is space-dependent. It is because 
in practice the quantum wells of different shapes are produced by means of MBE technique, when 
the semiconductor layers are grown subsequently one by one and these layers are characterized 
by different electron effective masses. The subsequent semiconductor layers stuck together make 
the quantum well and as a result, one can consider such structure as having effective mass which 
depends on the space variable. Proceed now to the thorough treatment of this problem. 




and thereupon, to a simple differentiation of its kernel 



V(x) = -2—K(x,x) 



Q(x, x') = — [1 - S(k)} exp(ik(x + x'))dk + ^ M% exp(-K n (x + x% 




n 
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3 Reconstruction of quantum well potentials using intertwining 
operator technique 



3.1 ISP and reconstruction of quantum well potentials 

The discussion presented above suggests the following strategy of QW potential reconstruction. At 
first, one has to search for the solution of our problem among the potentials of a special class of 
reflectionless and symmetric potential V(x) = V(—x), using the ISP-method as it was done in [4]. 
Then one can amend the potential reconstructed in this way, using the intertwining operator tech- 
nique and taking into account the dependence of the electron effective mass on the space variable. 
Now we are able to describe this recipe in more details. Suppose that one-dimensional potential 
V(x) can be represented by the function Vn (x, m*, Eq) which obeys the following conditions: 

(i) Vn supports precisely N bounded states of the quantum system with the effective mass m* . 
The bound-state energies coincide with the energies e±, £2, ...ejv of the levels within the QW ; 

(ii) Hindoo Vn = Eq. The last value can be considered as the depth of QW. 

Arranging the binding energies k 2 = Eq — e n in descending order so that k\ > > ••• > few, 
and e\ = Eq — k\ refers to the ground-state energy, one can use for QW potential reconstruction 
the technique developed by Schoncfeld et al [8] for studying the convergence of the reflectionless 
approximation to the confining potentials. Omitting the intermediate calculations, we give here 
only the final results: 



Here the sum ranges over all subsets S of {1, 2, N} including the null set and the full set, while 
S denotes the complement of the set S. 

Now we consider a generalized Schrodinger equation, that is with the potential Vn(x,Eq) and 
the position-dependent effective mass. Obviously, the new spectrum obtained in this way might 
differ from that which was used in ISP-method in order to reconstruct the potential Vn(x, Eq), but 
we can suppose that the changes of a spectrum are not dramatic, because in practice the space 
variable dependence of the effective mass is weak. Further on we shall show that it is possible to 
amend this potential and obtain "improved" one V(x), in order to have the spectrum needed. 

3.2 First-order Darboux transformation and supersymmetry 

The problem of the space-variable dependent effective mass attracts now persistent attention be- 
cause it is not obvious whether the effective mass approximation is applicable to heterostructures, 
or not (see [9,10] and the Refs. therein). Before we start our discussion, let us make some general 
remarks. Remember, the QW of the shape other than rectangular, is produced by stacking up a 
number (some times even a considerable number) of layers of different materials, each of which 
is characterized by its own effective mass. This stack of layers can be considered as some special 
case of heterostructure. Then if one tries to solve the Sturm-Liouville problem for a corresponding 
Schrodinger equation (we refer to this problem also as a direct one), treating the heterostructure 
as a whole and applying the effective-mass theory, one encounters some difficulties, whose nature 



V N (x,Eo) = Eo-2—~hiD(x) 



(2) 



where 




(-2xY,kp)l[(S,S) 




U(S,S) = 
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is the following. First, the total electron wave function is a product of the slowly varying envelope 
function and the Bloch function of the local extremum in the host's band structure. The Bloch 
functions in the two materials on either side of a heterojunction must be similar for the effective- 
mass approximation to be valid. An obvious condition is that they must belong to the same point 
in Brillouin zone, and this can fail for some materials. The second point concerns the matching 
of the envelope functions at the interface. Consider a junction at x = between two regions of 
materials, say A and B. The Schrodinger equations for the envelope function in the two subsequent 
regions (we consider only one-dimensional model), are 

1 d 2 



m* A dx 2 



m* B dx 2 



where m* A and m* B are the electron effective masses for the materials A and B, respectively, h 2 /2 = 1 
and the difference in the bottoms of the conduction bands is AE C = E B — E A . If the materials 
were the same, one can match the value and the derivative of the wave function at the interface by 
means of usual conditions 

* (0 - } = (0+) ' -*r\*=°- = 

where 0_ means the side of the interface in material A and so on. This simple condition is not 
correct for the heterostructure where the two effective masses are different, because it does not 
conserve current. A correct set of matching conditions is 

ch A (0 )-ch B (0 ) 1 d<t>A{ - X \ - 1 d4>B{x) \ 

The condition for matching the derivative now includes the effective mass. A more mathematical 
argument is that the matching condition which does not include effective masses assumes that the 
Schrodinger equation takes the form 

lr-M + V{x)ct>{x)=£<]>{x). 

m*(x) dx A 

This is not Hermitian (or Sturm-Liouville form) if m*(x) varies, and many of the crucial properties 
of the wave functions might disappear as a result. It seems that a consensus concerning the form 
of Hermitian kinetic energy operator in case of the position-dependent effective mass is already 
achieved among the specialists (see Refs. [9,10] and Sec. 4 below). In this section however, we 
are dealing with the intertwining relation technique, which is universal and apply this technique to 
construct a chain of exactly solvable Hamiltonians whose kinetic energy operators are not Hermi- 
tian, if m*(x) is space dependent. In the next section (Sec. 4), we also develop an approach which 
gives the possibility to treat the case of Hermitian kinetic energy operator on the same footing. 
Let us start with the equation: 

«*<«) -«««), H—-J i - ) * I + VM, (3) 

where m*(x) is a position dependent "effective mass" and V(x) is supposed to be equal to Vn(x, Eq). 
This equation is reduced to the generalized Schrodinger equation of the form: 

Tio<fro(x) = £m*(x)(f)o(x), Ho = —d 2 /dx 2 + v(x), 
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where v(x) = V(x)m*(x). In fact, it is the Schrodinger equation with linearly energy-dependent 
potentials. The Darboux transformations for the Schrodinger equations with variable values of 
energy and angular momentum were suggested in [11] and in a more general form in [12]. Then 
in Refs. [13,14] algebraic transformations have been elaborated for a Sturm-Liouville problem for 
studying phase-equivalent linearly energy-dependent potentials and for constructing exactly solv- 
able three-body models with two-central potentials. On the other hand, the intertwining operator 
method provides the universal approach to creating new exactly solvable models and can be ap- 
plied to the operators of a very general form ( see for example [15-17]). In this paper, we apply 
the intertwining operator technique to the equation (3) with a position-dependent mass in order 
to construct the potential which supports the desirable spectrum. 

Suppose that the solution of the eigenvalue problem to the equation ([3|) with the given potential 
V{x) and position dependent m*(x) are known and we would like to solve a similar problem for 
another Hamiltonian H containing a new potential V(x) and the spectrum which probably differs 
from the spectrum of the Hamiltonian © by a single quantum state: 



H<l>(x) = £<j){x) 



n 



i 



d 2 



+ V(x). 



m*{x) dx 2 

We start with standard intertwining relations (see, for instance [15,16]): 



(4) 



cn = hc, 



(j){x) = Of>(x), 



(5) 
(6) 



where the operator C intertwines the Hamiltonians TL and H. We search for the intertwining 
operator £ in a general form 

C = B(x)d/dx + A(x), (7) 

where A{x) and B(x) are to be determined. Once the operator C is known, the solutions <p can be 
obtained from (6) by applying C to the known solutions <j>. To find the explicit form of C, we use 
the equations ©,(4) and the intertwining relations (5), (6): 



1 



d 2 



m* (x) dx 2 
After some algebra we arrive at: 



+ V(x) C<t>{x) = C 



d 1 



m 



(x) dx 2 



+ V(x) (j)(x). 



- — (A"<f> + 2A'<f>' + A<j>") - —(B"(j)' + 2B >" + Bf) + V(A(j) + 
m* m* 



A 



.J-J' + vA+b(-± 

m J \ m* 

and finally, to the next system of equations: 



B 



(-)V 

\m* / 



+ v 4> + Vcj) 



— A + 2 — B 

m* m* 

1 / 1/, 

— 2A + — B 

m* m* 



v m / m 



VB 



-BV, 



—a" + va = av + bV 

m* 



From (8) it immediately follows that 

2B'/B = -m*'/m*,B = C/Vrn* 



(8) 

0) 
(10) 

(11) 
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where C is an arbitrary constant. Prom (9), (10) one gets 

1 B" 1 2A' , . 

V = V + — + — "5" 12 

m* B m* B 

and 

A" + ( — 2A' + B")B~ l A = BV' . 

m* m* 

In order to integrate the last equation, let us introduce a new auxiliary function K(x) defined as 
A{x) = B(x)K(x). Then we arrive at a nonlinear differential equation 



(-K" + 2K'K - V'm*) +^~(k 2 - K 



0. 



Taking into account the relation V = v/m* and the first of the relations (11), the last equation can 
be easily transformed into another one, in a single unknown K only: 



-K + 2K K - v ) - — ( —K + K' z - v ) = 0. 

m* 



This one can be rewritten as 



which means that 



d f 1 

dx \ m* 



-K' +K 2 -v) ) =0, 



(1/m*) [-K +K 
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where v is an integration constant. The last equation is analogous to Riccati equation. Introducing 
a new function U(x) as K = —IA IA~ 1 and changing v = —A, one arrives at the equation 

- (l/m*{x))U"{x) + V{x)U{x) = \U{x). (13) 

Here U(x) is supposed to be invertible at all x. The last equation then is nothing else but the initial 
equation © which is supposed to be solved and £ = X is a point of spectrum of H. Therefore, we 
assume that the solutions of (13) are known for the given values of A. Having found the explicit 
form of B (see (11)), using the formula for K mentioned above, from the relation A = BK one 
gets A(x) = —C (lnU(x)) \J\jm* (x). Once U is known, the transformation operator £, the new 
potential V{x) and the corresponding solutions of the transformed equation (4) are defined up to 
an arbitrary constant C. Without loss of generality, we can put it safely equal to unity. After this 
we have 



B(x) = l/Vm*R, A(x) = K/y/m*(x), K = - (\nU{x)) . (14) 

To make further transformations, let us calculate B" /B = \Jm*{\/ 'vra*)" '. Using this and (I14p in 
(7), (12) and (6) we construct the intertwining operator C, the transformed potential V(x) and the 
solutions (h in the form: 



(15) 



V = V + 



dx 1 



+ 2 



"rrr 



dx 



'rrr 



-.K 



V + 



"rrr 



dx 2 



"m* 



dx 



"rrv 



-.(Inn)' 
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Ccf> 1 



(17) 



"mr 

It follows immediately from (17) that CIA = 0. In order to obtain the solution of the equation (4) at 
the energy of transformation A, we shall use the second linear independent solution to (J3j) ^ namely 
lA(x) = U{x) f x dx'\U(x')\~ 2 where the integration limits depend on the boundary conditions. In 
particular, for regular solutions satisfying the boundary conditions <f)(x = 0) = 0, <p'{x)\x=o = 1 ; 
the lower integration limit is and the upper one is x, while for the Jost solutions the integration 
limits are — oo and x, respectively. As a result we get 

" w = ^ = jmm (18) 

Once 7] is found, one can get a second solution of (4) at the energy of transformation A. By using 
the Liouville's formula once more, one gets 

/x -i rx 
dx'\Tf\~ l = —== / dxlA{x)m*{x')lA{x'). (19) 
y / m*(x)U(x) J 

Hence, the information about all solutions of the initial equations ([3]) provides the knowledge of all 
solutions of the transformed equations ©. As in the case of Schrodinger equation, the functions 
<j)(x,£) and (p(x,£) correspond to Hamiltonians TL and Ti., respectively, are related through the 
transformation operator C (see (17)). The difference is that in our case C includes the position- 
dependent mass. As a consequence, the new potential V and solutions <f> depend on the effective 
mass m*(x). The function r](x) defined by (|18|) at the energy of transformation £ = A cannot be 
normalized and this is the reason why A does not belong to the discrete spectrum of 7i. Therefore, 
Hamiltonians H and H are isospectral with one exception of the bound state with the energy 
£ = A, which is removed from the initial spectrum of TC. Note that if the transformation function 
U(x) corresponds to the ground state, i.e., U(x) is nodeless, then the transformed potential V(x) 
has no any new singularity, exept the singularities due to V(x) (of course, we asumme m*{x) ^ 
at all x). However, if we apply this transformation to an arbitrary state other than ground 
state, the transformed potential V(x) might contain extra singularities, which are not present in 
the initial potential V(x) and hence, the Hamiltonian H becomes physically meaningless. As we 
shall see later, the difficulties with singularities can be circumvented by means of second-order 
Darboux transformations. Now we show how one can construct a Hamiltonian with an additional 
bounded state with respect to the initial Hamiltonian by using factorization of Hamiltonians and 
supersymmetry. 

The supersymmetry is based on factorization properties of Darboux transformation operators C 
and C + . The definition of formally conjugate operators is = (CQ)^ = Q^C^ and (^)^ = — J^- 
In our case, the scalar product of functions is defined by not the standard way (/, g) but with the 
weight of m*(x): (f,g) m = / m *( x )f{ x )9{ x )- I n this case instead of operator it is necessary 
to consider the operator m*~ l D^m* . Therefore the operator C) adjoint to C = -j={^ + K) is 
determined as 

£t * ( ' £ + A (20) 



"m 



dx 2m* 



Now let us consider the superposition DC, and CC^: 



£t £ = -—J- + —A-K' + K 2 ), (21) 
m* ax z m* 

Id 2 1 , . o . lm*" 3m*'m*' m*' , , 

££t = + - :{K' + K 2 - -— ^ + - -= -zK. 22 

m* ax z m* 2 m* z 4 m* J m* z 
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Express the potential V from equation ([13]) in the form V = U" / (m*U) + \. Using K' = —[U'/U]' = 
-U"/U + (U'/U) 2 we represent V as 

V = — (-K' + K 2 ) + \ . (23) 
m* 

Substitution of (|23j) into (116[) leads to the following representation of the transformed potential: 

V = —^K' + K 2 ) + - ^-K + A . (24) 

m* ym* V m * m 

Using (|23p and (|24p , after some transformations the formulae (j2T1 and (|22p can be rewritten as 

^C = ~^ + V-X = H-X; (25) 
1 

££ f = t^t + U- A = ft-A. (26) 

m* ax A 

/,From (|26p one can obtain the intertwining relation 

HC ] = C ] H, (27) 

which means that the operator is also the transformation operator and realizes the transfor- 
mation of the solutions of equation dH) to solutions of ([3]), cf) oc JUcfi. As one can see from the 
comparison of the relations (|15p and (|20p . the operator is not an inverse of C. One can show 
that the operators C and C) can be expressed in terms of 77, which are solutions of transformed 
equations (jU) at the energy A with the potential V determined by (fT6|) . For this aim let us express 
K in terms of rj, by means of (|18p . 

_ W _ m*' 7]' 
U 2m* rj 

Using this in (|15p and (|20p . we obtain 



£ _^f- + iSL + i), £ t_ _L f-^+a!] ( 2 8) 

ax 2m* 7] I ym* \ dx rj 



"m 



Evidently, the function 77 is also a transformation function. It is clear that £)r\ = 0, i.e., rj belongs to 
the kernel of the operator C) . As one can see from (|28p and (I19p . the application of the operator C) 
to the second linearly independent solution 77 to equation ^ gives back the solutions U of the initial 
problem at the energy of transformation. Indeed, C)f] = ^— ^ + ^\ i](x) j x dx'lrj 2 ^ 1 = U. 
Hence, a one-to-one correspondence between the spaces of solutions of equations ([3]) and (JU) is 
established, and these are the operators £ and C\ which produce the correspondence. 

Note, one can interchange the role of the initial and final equations. The function 77 becomes 
transformation function for the intertwining operator C\ which will make the transformation in 
the opposite direction: from the potential V to the potential V and from the solutions of to the 
solutions of ([3|). So, if within the first procedure {[T5j)— (JTTJ) we constructed the potential V with one 
bounded state removed , now we can construct the potential V with an addional bounded state. 

3.3 Second-order and the chain of Darboux transformations 

Let us define the second-order Darboux transformation as a sequence of two Darboux transforma- 
tions performed in a row 

C = C 2 £i, (29) 
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where C\ is actually C defined in (14) 



l 



d 

- r +K 1 

m* \ dx 



whereas Co is determined as follows: 



£ 5 



"m 



dx 



Ui' 



A 

X i 



(30) 



(31) 



and xi = Xi( x j^2) is obtained by means of the first-order transformation, applied to the solution 
U 2 of the equation (13) or © with the eigenvalue A2 



Xi = CiU 2 



1 (LJL w , 



"m 



* V dx hi 1 



(32) 



It is clear that xi is the solution of equation (13) with the potential V\ = V + 2K[, defined as in 
(|16p . and xi can be taken as a new transformation function for the Hamiltonian 7i\ to generate a 
new potential 



V 2 = Vi + 



and corresponding solutions 



y/m* (x) 



£■291 



d 2 



1 



d 

+ 2— 



1 



dx 2 yjm*(x) dx \ y/ m *(x) 



1 



d „ 
h K 2 

rrr \ dx 



>1, 91 



Cic 



(33) 



(34) 



Here the function (f>i, denoted earlier as 0, is an eigenfunction of the Hamiltonian 7i\ 

d 



>m 



(35) 



In other words, the action of the second-order operator (|29|) on the solutions <p leads to the solutions 

of n 2 

fa = C4 = C 2 C lV . (36) 

Iterating this procedure m times in regard to given operator Tt, one arrives at the operator H. m , 
which satisfies the intertwining relation 

CH = H m C. 



In this way one gets 



— ^m— 1 "I" 



1 



"m 



d 2 1 d 
= + 2— 



dx 2 ' dx \ v m 



1 



m—1 



9m — £>4> — £-m9m-l — C, m L m -\...L\l. 

where C is the m-th order differential operator: 



(37) 
(38) 

(39) 



It should be noted, that the chain of m first-order Darboux transformations results in a chain of 
exactly solvable Hamiltonians Ti — * TL\ TC m . 
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Consider now the 2-nd order transformation in detail. Using the explicit expression for V\ which 
appears in the first-order Darboux transformation, we get a formula for the potential V 2 : 



V 2 = V + 



"m 



d 2 



1 



dx 2 



"m 



+ 



d 



'm 



* dx 



1 



"m 



-.K 



where K = K\ + K 2 . Let us represent Xl as 

Xx(x) = 



W 1>2 (x) 



(40) 



(41) 



■sjm*(x) Ui(x) 

where W\ t2 {x) = U\{x)U' 2 [x) — U[(x)U 2 (x) is the Wronskian of the functions U\(x) and U 2 (x). 
Plugging (|4ip into the formula (|3ip for K 2 , after some transformations we obtain 



m 



W 1)2 (x) 



yJm*{x)Ui(x) 



(42) 



After this K 



K\ + K 2 can be represented as 



K 



U[ m*_ 
Tlx + 2m* 



+ 



in 



U x W 1)2 2m* W 1>2 
With this taking into account, making in (|40p the next substitution: 

ldm' d 2 1 



2 dx m* 3 / 2 dx 2 m* 1 / 2 ' 
after some manipulations the new potential can be expressed as: 



Vo(x) 



V(x) 



d 



dx 



1 d 



'm 



dx 



InWi, 2 (x) 



(43) 



By using (f34l) find now the corresponding functions <j) 2 (x), <f> 2 = + K 2 ) By analogy with xi 
the function <t>\(x) can be written in terms of the Wronskian Wie(x) = U\{x)<$ {£ , x)— U[(x)(f)(£, x): 

1 Wx, £ (x) 



cj>i(x) 

Let us now calculate the derivative of 



y/m*(x) U\(x) 
C\4>, that is 

1 ^// 1 W i 



(44) 



+ 



/ m*W 1 \Jm* \Jm* U\ 

Making use of the last expression and the relation (|42p for K 2 , we obtain, after some simplification, 
the formula 



(j) 2 (x) 



1 



m*(x) 



fix) 



m*(a;)Wi(a;) 



(45) 



It is easily seen from (|43p and (|45|) that due to the 2-nd order Darboux transformation, the potential 
and solutions obtained in this way are completely expressed in terms of the known effective mass 
function m*{x) and the solutions U±(x) ,t( 2 (x) , 4>(£, x) to the initial equation, with no use of the 
solutions to the intermediate one with the potential V\{x). 

Clearly, for the next transformation step to be made, one should take a new transformation 
function X2, that corresponds to the potential V 2 . It can be obtained by applying the operator 
C = £2^1 to the solutions U3 corresponding to the eigenvalue V3: 



1 



X2 



m*(x) 



^(lnW 1 , 2 (a)) 



W 1>3 (x) 



m* (x)Ui(x) 

Then it can be used to produce a new transformed operator £3 = d/dr + K3, 
generating new potential V3 and solutions cps and so on, according to (|37|) -(p9 



~X 2 X 2 



for 
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3.4 The integral form of Darboux transformations 

The transformed solutions (|44p and (|45p can be represented in the integral form. Let us consider 
to this end the generalized Schrodinger equation written down as 

- (j)"{x) + m*{x)V{x) = £m*(x)<t>(x). (46) 

Multiplying both sides of the equation (f46l) for the function 4>(£,x) by U\(x) at the energy of 
transformation Ai and subtracting from the obtained expression the equation similar to (|46|) but 
written down for li\{x) and multiplied by (j)(£,x), we arrive at 

■fwiA*) = ( A i " £)m*(x)Ui(x)<f>(£,x). (47) 
ax 

The last expression can be easily integrated: 

W he (x) = (Ai - £ ) f m*{x')Ui{x')^{x')dx' + C. (48) 



Inserting the last expression into the formula for <j)\ (|44p . we arrive at the integral form of the 1st 
order transformed solutions: 

[C + (Ai - g) m^x'^jx'^jx^dx'} 
Mx) = b f ; ■ (49) 

Here C and a are some arbitrary constants. By analogy, applying this technique to the equation (J3j) 
for <fi and Ui, using (|48p in (|45p. we get the integral form for the 2-nd order transformed solutions 

, a c W ^ d fl W ( ^^ ( g + (Al ~ £) !l rrfVnWtWM) 
H = (A X -mx) --{toWwix)) m * {x)Ulix ) • ( 5 °) 

Here the integration limits depend on the boundary conditions. In particular, for regular solutions 
satisfying the boundary conditions cj>(x = 0) = 0, (ft'(x)\ x= Q = 1, the lower integration limit is and 
the upper one is x, while for the Jost solutions the integration limits are x and oo, respectively. 
The constant C is determined by the values of the Wronskian at zero or at infinity, depending on 
the way the problem is posed. Notice that the functions U and (f> can be chosen in such a way that 
the constant C becomes zero. Analogously to (@7J), one has W[ 2 (x)/(\i — A 2 ) = m* (x)Ui(x)U2{x) 
and 

W 1>2 (x) ci + j x dx'm^x^U^x'^ix')' 1 ' 

Using the last formula and assuming C = 0, after some transformations one can represent (j>2 as 
follows: 

<p 2 = (Ai - £)4>{x) TTWrf Vr a77 7 a77 7 n ■ ( 52 ) 

ci + J ax'm*{x')Ui(x)U2{x) 

Now let us consider the 2-nd order Darboux transformation at Ai = A2 = A. Earlier within 
the first-order procedure, we already obtained two linear independent solutions (|18p and (I19p at 
Ai = A2. The second transformation can be made by means of a linear combination of the solutions 
rj and fj 

Xi(x) = cir](x) + rj(x) = / 1 „ </ — ■ ( ci + / dx'U 2 (x')m*(x') ) . (53) 

^m*(x)U(x) V J J 
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In order to find the transformed potential and solutions, calculate K2 = — Xi/Xi an d K = K\ + K2 

m*' (x) m*(x)U 2 (x) 



K x) 



2m* (x) (ci + f x dx'U 2 (x')m*(x')) ' 
Plugging the last expression into the formula (I40|) which defines the potential, we arrive at 



M M 2 d 1 ^ 2 (x)m*(x) \ 

^ " [X) ^{x)dx \^ m T(^(c 1 +J x dx'U 2 (x')m*(x') J' 1 j 



The operator £2 flM]) with xi defined by (|53|) . acts on the function </>i represented by its integral 
form (|49p so that it leads to 

ch = (a - - ^)( A -£)r^wv*w(^ (55) 

ci + J* dx'U 2 (x')m*(x') 

It is worth mentioning, that the formulae for the new potential V2 and the solution <p2 can be 
obtained directly from the relations (03]) and (f4"5|) . if one takes into account that at Ai = A2 = A, 
the expression (j5Tj) for ^ In 1^1,2(2:) should be changed by 

d .lnP(x)- -W 2 (*) 



cix ci + / x dx'U 2 (x')m* (x') ' 

with P(a?) = ci + J x dx'U 2 (x')m*(x'). 

Without loss of generality one can take the linear combination of the functions 7/ and fj as 
Xi(x) = 7/(x) + Crj(x), and change (A — £)4>{x) — ► for simplification. Then formulae ([54"]) and 
(|55p can be rewritten as 

yn _ yn 2 d / 1 CU 2 {x)m*{x) \ 

2{X) [X) ^{x)dx \^^)(l + CjUx'W{x')m*{x') ) ■ [M) 

U{x)Cf x dx'U(x , )m*(x')<K X ') 
^ = ^~ l + C f: o dx>U 2 (x>)m*(x>) ■ (57) 

The constant C plays the role of a normalization constant or the difference between the normal- 
ization constants of the bound state A for the potentials and V(x), respectively. Notice, the 
choice of arbitrary constants x Q and C allows one to avoid the problems with zero-equal denom- 
inators, or in other words, it means that one can make transformations on an arbitrary bounded 
state and construct the potential without singularities. Notice also, that m*(x) itself does not lead 
to the singularities, because the effective mass m* (x) / and assumed to be smooth and at least 
twice differentiable function with respect to space- variable. 

The solution of the equation ([3]) with the potential (|56p at the energy of transformation A can 
be achieved by means of operator £2 acting on the solution 77 from (|18p . obtained within the first 
transformation step 

r]2{x) = C 2 r) ~ 



s/m*{x) \dx x{x) J ^m*(x)U(x) 

where x' is assumed to be of the form (|53p . Finally we get 

, s CU(x) 

V2[ - X> 1 + C j x dx'm*(x')U 2 (x'Y [ ' 
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One can rewrite the potential (|56p and the solutions (|57p in terms of r]2(x) as 



V 2 (x) = V(x) + f l_ ^y m *{x) m (x)U(x)] , (59) 

/x 
dx'U(x')m* (x')4>(x'). (60) 

The relations (|56p - (|60p are the results of performing two subsequent transformations with the 
same energy. Therefore, it allows one to construct the phase-equivalent potentials. Indeed, if 
C = N$ — N 2 is the difference in normalization constants of the bound state A for the potentials 
V2(x) and V(x) respectively, then the formulae ([56]) . ([57]) and ([58]) correspond to phase-equivalent 
potentials whose scattering data coincide and differ only by a normalization factor. Note, the 
phase-equivalent potentials have a different shape. They can be more deeper and narrow or more 
shallow and wider and possess the same spectral data, except for normalization constants. 

If we assume the transformation function U{x) to be taken at the energy of the bounded state, 
which we would like to add to the initial spectrum, and C = N 2 is the corresponding normalization 
constant, then the formulae ([56]) . ([57]) and ([58]) give the possibility to construct a potential with 
a new bounded state A provided the other spectral characteristics of the spectra produced by the 
potentials V^(x) and V(x), coincide. Notice, that the function U(x), which is the solution of the 
initial equation with the potential V, has to be taken at the energy of transformation A. To 
sum up, it can be said that by means of the technique described above, it is possible to remove 
some bounded y states or to add new ones and to construct the phase-equivalent potentials. The 
procedure can be repeated as many times as it is needed to construct a new potential with a 
desirable spectrum. Finally, it should be noted that using the technique presented in [12,13], it is 
not difficult to generalize these results to include a case when the second order transformation is 
applied simultaneously to N bounded states. This, however, is beyond the scope of the paper and 
will be a subject of another publication. 



4 Liouville transformation and QW potential reconstruction 

In previous sections we have reconstructed the QW potential in the downright fashion using, in fact, 
the position-dependent potential function. That is why we did not worry about the hermiticity of 
the kinetic energy operator. 

In this Section we propose a somewhat different approach to the reconstruction of QWs poten- 
tials, which is based on the following observations. First, in most of cases the number of bounded 
states in QW is limited from above by some reasonable value of about 6-10 which, in its turn, 
is determined by the technological limitations. Secondly, in most of cases one can safely suppose 
that effective mass m*(x) is a smooth and at least twice differentiable, weakly position-dependent 
function. Therefore, we can use the Liouville transformation together with algebraic Darboux 
transformations in order to reconstruct the QW potentials with a predetermined spectrum. This 
approach is applicable to the case of the straightforward reconstruction of QW potential discussed 
previously, as well as to the case when the Hermitian kinetic energy operator enters explicitly the 
generalized Schrodinger equation. 

Let start from the equation of the form 

--^<%) = m*(g)A<%), (61) 
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where q stands for space variable. Thus, one can interpret this equation as a generalized Schrodinger 
equation with the potential V = 0. Introduce now new variables and a new function as follows 
(this is one of the special cases of Liouville transformation ): 



m = p , \mr 



dx 
dq 



, $( g ) =p- l / 2 m{ x ). 



(62) 



Then, taking into account ([62]) . the equation (f6"Tj) can be reduced to the following: 



-P 



3/2 



* (x) + 



1)1 E* 

2 | 2 I p 



p 2 Xp 1/2 ^(x 



which can be written down as a standard Schrodinger equation 

-^"(x) + V(x)^(x) = A*(x), 
where the potential V(x) is of the form: 



V(x) 



1 



Px_llP JL 
P 2 I p 



, \ 2 



Since 



p dx \ p J \ p 



/\ 2 



the potential V(x) can also be represented as 

1 



^(x) 



or returning back to the variable m* , 



d 2 1 1 / d . 
dx^ lnp+ 2W nP 



(63) 



V(a;) 



m 



3 m* 

4 \ m* 



Now assuming the solution of the eigenvalue problem to the Schrodinger equation with this potential 
is known, one can add or subtract one by one as many bounded states as it is needed, using the 
algebraic Darboux transformations described in the previous Section. 

As for the effective mass theory with the Hermitian kinetic energy operator is considered, it 
seems that a consensus among the specialists has been already achieved, since in the majority of 
corresponding studies the following form of the Hamiltonian is used (see, for example, [10]): 



H 



P 



1 



m*{q) 



P 



+ V(q), 



where P is the momentum operator. Then, in order to use the Liouville transformation approach, 
one should start with the equation 



d 



1 d 



dq \ m* (q) dq 



15 



Making the next substitutions 

— 1 _ dq _ _i/2^ 
y/m*(q) dx' P 

after some manipulations similar to that which were previously done, one get the standard Schrodinger 
equation 

-^"(x) + V(x)W(x) = A*(x), 
with the potential defined as in (f63|) but with p{x) = m*~ 1 (x). Finally we obtain 



V(x) = \ 



m*" 5 / m*' \ 
m* 4 \ m* I 



5 Conclusion 

The basic elements of contemporary micro- and nanoelectronics are the low-dimensional structures 
which are the structures composed of QWs, quantum wires and quantum dots and produced by 
means of various techniques including most impressive one, molecular beam epitaxy. The entirety of 
such methods and techniques are sometimes termed as Quantum Engineering or Quantum Technol- 
ogy. One of the most important issues of quantum engineering is the construction of multi-quantum 
well structures possessing desirable properties. This problem appears in different contexts, ranging 
from the construction of multi-level computer logic to photovoltaics of third generation [18,19]. 
From the theorist's point of view, the problem can be formulated as follows: assume one requires a 
definite spectrum of QW, because it is determined by some specific needs and circumstances. Can 
one reconstruct the QW potential which supports this very spectrum? In this paper we answer 
this question in affirmative and outline the possible strategy of the QW potential reconstruction, 
if the spectrum of QW is predetermined. 

The proposed approach is based on the combination of various techniques such as Inverse 
Scattering Problem Method, Darboux and Liouville transformation. Bearing in mind that the ef- 
fective masses of charge carriers in the subsequent layers of different materials which make QW, 
are different, we match the intertwining operator technique, in order to take into account the 
position-dependent mass in Eq. (j3J). The first- and second-order of Darboux transformations, as 
well as the chain of Darboux transformations are considered, and interrelation between the differ- 
ential and integral transformations is established. The developed approach allows one to construct 
phase-equivalent potentials and to add (or if necessary, to remove) some states to (or from) the 
spectrum supported by the initial potential, whose form can be established for instance, by means 
of ISP-method. 

Another possible way to take into account the position-dependent masses of charge carriers is 
to use the Liouville transformation. It also allows one to incorporate and treat on the same footing 
the proper form of Hermitian kinetic energy operator. This operator appears in the context of 
applicability to heterostructures the effective mass approximation, the subject of today's research 
activities. 

In the paper we only formulate the basic concepts and develop the necessary mathematical 
tools. More thorough numerical examinations of specific cases are required and it will be done 
elsewhere. 
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